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ABSTRACT 

We use three dimensional magnetohydro dynamic (MHD) simulations to model the 
supernova remnant SN 1006. From our numerical results, we have carried out a polar¬ 
ization study, obtaining synthetic maps of the polarized intensity, the Stokes parameter 
Q, and the polar-referenced angle, which can be compared with observational results. 
Synthetic maps were computed considering two possible particle acceleration mecha¬ 
nisms: quasi-parallel and quasi-perpendicular. The comparison of synthetic maps of the 
Stokes parameter Q maps with observations proves to be a valuable tool to discern 
unambiguously which mechanism is taking place in the remnant of SN 1006, giving 
strong support to the quasi-parallel model. 

Key words: MHD-radiation mechanisms: general - methods : numerical - super¬ 
novae: individual: SN 1006 -ISM: supernova remnants 


1 INTRODUCTION 


Recently, the study of bilateral (also known as barrel-like) 
supernova remnants (SNR) has gained great interest since 
they have proven to be a useful tool when studying the con¬ 
figuration of the interstellar magnetic field (ISMF) on scales 
of a few pc. As their name suggests, these type of SNRs dis¬ 
play two characteristic bright and opposite arcs. The mor¬ 
phology of this type of remnants has been explained with 
an ISMF, which is involved in the acceleration mechanism 
of relativistic particles. However, there is still debate as to 
which is the most efficient acceleration mechanism in these 
objects. One interpretation is the equatorial belt model, in 
which the orientation of the ISMF is quasi-perpendicular to 
the shock front normal. The other interpretation is the polar 
cap model in which they are quasi-parallel. 

The remnant of SN 1006 is the archetype of the bilateral 
SNR group. It has a diameter of 30'(or 20 pc at a distance 
of 2.2 kpc, Winkler et al. 2003) and, in radio-continuum 


and X-ray images, exhibits an incomplete shell with two 


bright arcs perpendicular to the Galactic plane (Reynolds 
&, Gilmore||1993| Reynoso et al.|[20l3 ). SN 1006 is accepted 
to be the result of a type la SN explosion (Stephenson & 
Green|2002 ). This remnant is located at a high Galactic lat- 
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itude, and therefore thought to be evolving in an almost ho¬ 
mogeneous interstellar medium (ISM). For this reason, the 
radio-continuum morphology of the remnant must be mostly 
affected by the characteristics of the surrounding ISMF. 


Numerous theoretical and observational studies have 
been devoted to analyze which mechanism is responsible for 
the morphology of this remnant, giving rise to two opposing 
interpretations: either the quasi-perpendicular acceleration 


mechanism (e.g. Fulbright V Reynolds 1990 Petruk et al. 

[2009 

Schneiter et al. 2010) or the quasi-parallel mechanism 

(e-g- 

Volk et al. 2003 Rothenflug et al. 2004 Bocchino et al. 

2011 

)• 


Most of these works apply criteria based only on the 
brightness distribution of the synchrotron emission. Accord¬ 
ingly, their conclusions are strongly dependent on the ori¬ 
entation of the ISMF respect to the fine of sight ( los ), and 
can not determine which particle acceleration mechanism 
is actually taking place. A recent observational polariza¬ 
tion study seems to favour the idea that the bright arcs 
of SN 1006 can be explained by the polar cap model, imply¬ 


ing an ISMF parallel to the Galactic plane (Reynoso et al. 


2013). 


The present work is an effort to clarify which process or 
mechanism is more suitable to explain the synchrotron emis¬ 
sion of the remnant of SN 1006. For this purpose, we carried 
out 3D MHD simulations employing the same scenario pro- 
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Figure 1 . Scheme of the numerical setup employed in our 
simulations for the quasi-parallel (top panel) and the quasi¬ 
perpendicular case (bottom panel). In all the models a dense cloud 
was placed to the right. 


posed in Schneiter et al. (2010). From the numerical results, 
we performed a polarization analysis of the Stokes param¬ 
eters Q and U, as well as the polar-referenced angle, thus 
facilitating the direct comparison with observations. 

The organization of this work is as follows: section 
presents the initial setup for the MHD simulations, section 
^3] explains how the synthetic maps for the radio emission 
and Stokes Q and U parameters were calculated, section © 
introduces the results obtained which are further discussed 
and summarized in section ^51 


The SNR 


The physical parameters of the SNR setup are the same as 
those presented in Schneiter et al. (2010), but to help the 


readability of the present paper, we reproduce them below. 

A supernovae explosion is initialized by the deposition 
of E 0 = 2.05 x 10 51 erg of energy in a radius of Ro = 0.65 
pc located at the centre of the computational domain. The 
energy is distributed such that 95% of it is kinetic and the 
remaining 5% is thermal. 

The ejected mass was distributed in two parts: an inner 
homogeneous sphere of radius r c containing 4/7ths of the 
total mass (M* = 1.4Mq) with a density p c , and an outer 
shell containing the remaining 3 /7ths of the mass foll owing 
a power law (p oc r ~ 7 ) as in Jun h Norman (1996b). The 


velocity has a linear profile with r, which reaches a value of 
vq at r = Ro. The parameters p c , r c , and vo are functions 
of Eq , M*, and Ro, and were computed using Eqs.(l)-(3) of 


Jun & Norman (1996b) 


The surrounding interstellar medium 


The surrounding ISM, where the SNR evolves, was simu¬ 
lated as an homogeneous plasma with temperature To = 
10 4 K and number density no = 5 x 10 -2 cm -3 . At a dis¬ 
tance of ~ 8pc in the ^/-direction from the centre of the 
computational domain, the density was increased by a fac¬ 
tor of 3 ( Acero et al.|2007 ). The SNR will collide with this 
‘wall’, producing a bright Ho emission filament (Winkler 


et al. 

2003 

, which has X-ray and 

et al. 

2007 

Schneiter et al.|2010) 


In order to simulate the synchroton emission of SN 1006, 
based on the results of Schneiter et al. (2010) we have con¬ 


sidered two configurations for the magnetic field B depend¬ 
ing on which of the particle acceleration processes is as¬ 
sumed. For the quasi-perpendicular case, B is along the 
x — direction, and it is along the y — direction for the quasi¬ 
parallel case (the bright arcs are in the y — direction). In 
both cases, we have considered a mean magnetic field mag¬ 
nitude of 2/xG, while only for the quasi-parallel case an ad¬ 
ditional VB, increasing toward the —x — direction, was in¬ 


cluded (which corresponds to model GRAD2 of Bocchino 
et al.|2()T l ). Figure [l] shows an scheme of the initial numer¬ 
ical setup employed in our simulations. 


SYNTHETIC TOTAL AND POLARIZED 
RADIO EMISSION MAPS 


2 THE NUMERICAL MODEL 


2.1 Initial setup 


The numerical simulations were carried out with the paral¬ 
lelized 3D MHD code Mezcal ( De Colle &; Raga||2006| De 
Colle et al.|2008 2012). This code solves the full set of ideal 
MHD equations in cartesian geometry (x',y',z') with an 
adaptive mesh, and includes a cooling function to account 
for radiative losses (De Colle & Raga 2006). The computa¬ 
tional domain is a cube of 24 pc per side and it is discretized 
on a six level binary grid, with a maximum resolution of 
2.3 X 10 _ 2 pc. All the outer boundaries were set to an out¬ 
flow condition (gradient free). 


3.1 Synchrotron emission maps 


The synchrotron emission is obtained as (Ginzburg &; Sy- 


rov atskii|' l965): 


i{y) oc KB° ± 


( 1 ) 


where v is the observed frequency, B± is the component of 
the magnetic field perpendicular to the los , and a is the 
spectral index which was set to 0.6 for this object. The coef¬ 
ficient K oc e v~ h includes the dependence from the obliquity 
and the shock velocity v 3 (O rlando et al.| |2007), where e is 
either proportional to sin 2 O b s for the quasi-perpendicular 
case or cos 2 O b s for the quasi-parallel case. The angle O b s 
is the angle between the shock normal and the post-shock 
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Figure 2. Synthetic maps of polarized intensity obtained for the 
quasi-parallel (panel b) and quasi-perpendicular (panel c) cases. 
The axes are in pc and the linear color scale is in arbitrary units. 
These maps were rotated 60 degrees with the purpose of com¬ 
paring them to a radio image at 1.4 GHz (panel a), displayed in 
equatorial coordinates. 


magnetic field ( Fulbright &; Reynolds! 1990 ). As in _ 

et al. (2007), the exponent b was chosen to be —1.5, implying 


Orlando 


that stronger shocks are more efficient injecting particles. 


3.2 Maps of the Stokes parameters 


To calculate the Q and U synthetic maps, the following 
expressions for the Stokes parameters were employed (see 
Clarke et al.||1989 Jun &; Norman) 1996a ) 



foi(v) cos [2<j>(s)] ds 


foi(v) sin [ 2 4 >(s)] ds 


( 2 ) 

(3) 


where s increases along the los , 0(s) is the position angle 
of the local magnetic field on the plane of the sky, and /o 
is the degree of linear polarization, which is function of the 
spectral index a: 


n _ 0+1 

/o ” o + 5/3 


(4) 


The linearly polarized intensity is given by: 

MO = RQM 2 + MM (5) 

and the polarization angle was computed as follows: 
X=ltan ~\U/Q) (6) 


3.3 Observations 


In order to compare our simulations with SN 1006 as ob¬ 
served in radio wavelengths, we have made use of data ob¬ 
tained at 1.4 GHz with the Australia Telescope Compact 
Array (ATCA) and the Very Large Array (VLA), already 
published in Reynoso et al. (2013). These data contain in¬ 


formation on the polarization of the emission, thus we con¬ 
structed a map of the polarized intensity (Fig. If) and an¬ 
other map with the Q parameter (Fig. If) . The observa¬ 
tions and data reduction process are described in [Reynoso] 
et al. (2013). The images presented here are convolved with 
a beam of 15 arcsec. To remove spurious features beyond 
the limits of SN 1006, only pixels where the total intensity 
is above 75 mJy beam -1 are retained. 


4 RESULTS 

4.1 Synthetic polarization maps 

Synthetic polarized intensity maps were constructed from 
our numerical simulation results, using Eq.|5]). Figure[2]dis- 
plays a comparison of the observed polarized emission (panel 
a), and the synthetic maps for the quasi-parallel (or polar 
cap model, panel b) and quasi-perpendicular (or equatorial 
belt model, panel c) cases. Before integrating along the los 
the computational domain (the x y z' — system) was tilted 
60° counterclockwise around the z—axis (the /os), such that 
the North points upwards and the East to the left (the 
plane of the sky is the xy plane). In addition, both mod¬ 
els were rotated 15° around the y— axis. The quasi-parallel 
model was further rotated 15° in the a:—axis. Both synthetic 
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Figure 3. Comparison between the synthetic Stokes parameters 
Q for the quasi-parallel and the quasi-perpendicular cases (panels 
(b) and (c), respectively) with observations (panel a). 


maps exhibit two bright arcs in the NE and SW direc¬ 
tion, resembling the observations. The magnetic field gra¬ 
dient along the x —direction introduced in the simulations 
produces a SW-NE asymmetry in the polar cap model, in 
agreement with Bocchino et al. { 2011 } and with observa¬ 
tions (see Reynoso et al.||2013 and Figure [5] in this paper). 
As expected, the equatorial belt model shows greater inten¬ 
sity due to a higher efficiency in diffusive shock acceleration 


( Jokipii||1987 Rothenflug et al.||2004 Cassam-Chenai et al, 

20081 ). 


4.2 Stokes parameter Q 


In order to compare the simulated maps of the Stokes pa¬ 
rameter Q with the observational one, we had to take into 
account that the observation of radio polarization from the 
SNR actually gives the position angle of the electric field. 
Hence, we had to follow the inverse process to the observa¬ 
tional one, that is, the magnetic field was rotated clockwise 
90°, emulating the perpendicular electric field, and then we 
performed the inverse Faraday correction. To clarify, we have 
replaced the argument (f)(s) in the trigonometric functions 
of Eqs. © and © by (j)(s ) — 7t/2 + Axf, where A\f is the 
Faraday rotation correction, given by: 

* RM f A \ 2 

AXF “ [rad m- 2 ] ( [m] ) ^ 


being RM the rotation measure and A the observed wave¬ 
length. For the case of the SN 1006, Reynoso et al.| ([2013) 
obtained an average RM of 12 rad m~ z from their polariza¬ 
tion study at A = 0.21m. With these values, Eq. 0 gives 
a correction angle due to Faraday rotation of 0.52 rad or 
30.3°. 


Figure [3] shows the resulting synthetic maps of the Q 
parameter (panels (b) and (c) for the quasi-parallel and the 
quasi-perpendicular cases, respectively), which were com¬ 
pared to the observed map (panel a). 


As mentioned above, the question about which process 
of acceleration of relativistic particles is responsible for the 


synchrotron emission in SNRs (see |Fulbright <V Reynolds 

1990 

Petruk et al. 2009; Schneiter et al. 2010; Volk et al. 

2003 

Rothenflug et al.|2004||Bocchino et al.|2011|) is an on- 

going 

discussion in the literature. Our simulated polarized 


intensity maps are not conclusive on which of the two com¬ 
peting processes better explains the case of SN 1006. How¬ 
ever, the maps of the Stokes parameter Q are quite differ¬ 
ent between the two models, thus removing the discrepancy. 
Hence, the Stokes parameters turn out to be an important 
tool to compare numerical models with observations and de¬ 
termine which of the acceleration processes is more suitable 
in each scenario. 


4.3 Polar-referenced angles 


Following Reynoso et al. (2013), we created maps of the 
angular difference for the position angle x with respect to 
the local radial direction f, i.e. we obtained maps with the 
distribution of the polar-referenced angle % r , which is given 
by: 


Xr 




( 8 ) 
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where b± is the direction of the magnetic field on the plane of 
the sky. In this representation, Xr = 0° and Xr — ±90° cor¬ 
respond to the radial and tangential directions respectively. 
These maps allow us to analyze two aspects of the magnetic 
field: the orientation of the post-shocked magnetic field and 
the orientation of the unshocked magnetic held of the am¬ 
bient ISM. We expect the magnetic fields in young SNRs to 
be radially aligned due to the development of hydrostatic in¬ 
stabilities (Jun & NormanJl996a). So, the polar-referenced 
angle, Xr, should have a value close to zero on the edge of 
the shell. 

Figure [4] shows maps of the distribution of Xr, obtained 
for both quasi-parallel and quasi-perpendicular cases. These 
maps show a very small region on the SNR periphery with 
radial magnetic held orientation, and Xr takes values that 
run from 0°(blue color), to 90°(red color). For both cases, 
the regions with Xr — 0 are nearly aligned with the direc¬ 
tion of the unperturbed magnetic held as imposed in each 
model. This ‘radial’ region is tilted clockwise ~ 30° respect 
to the y or ‘North’ direction (see panels (b) of Figure |4| for 
the quasi-perpendicular case, and ~ 60° counterclockwise 
for the quasi-parallel case (see panel (a) of Figure [4|. The 
quasi-parallel model gives a very good agreement (at large 
scales) with the observations (see Figures 7(a) and (b) in 


Reynoso et al. 2013), suggesting that the direction of the 


ambient magnetic held proposed in this model must be rep¬ 
resentative of the actual orientation. We must note however, 
that the polar-referenced angles show some spread at small 
scales. This is most likely due to inhomogeneities of the mag¬ 
netic held at such scales, which is not accounted for in our 
simulations. 


5 DISCUSSION AND CONCLUSIONS 


We have utilized 3D MHD models to simulate the polarized 
emission of SN 1006. The 3D simulations allow us to explore 
configurations that were impossible with the axisymmetric 


code used in Schneiter et al. (2010). In particular it is pos¬ 


sible to test different orientations of the magnetic held with 
respect to the plane of the sky and avoids some artifacts 
that arise when generating synthetic maps. 

SN 1006 is a member of a population of ‘barrel-like’ 


or ‘bilateral’ SNRs. Previous studies (Reynolds &; Gilmore 


1993 IPetruk et al.|2009||Schneiter et al.|2010||Bocchino et al. 


2011 Re ynoso et al.|2013 ) have shown that the emission ge¬ 

ometry of bilateral SNRs is correlated with the orientation 
of the ambient interstellar magnetic held. Different orienta¬ 
tions of the magnetic held can be inferred depending upon 
which particle acceleration mechanism is assumed to be tak¬ 
ing place at the SNR shock front, either the quasi-parallel 
or the quasi-perpendicular model. 

As mentioned above, there is an ongoing debate as to 
which of these two particle acceleration mechanisms is re¬ 
sponsible for the observed syncrotron emission in SNRs. 
While observational works seem to support the quasi¬ 
parallel model (polar cap scenario), theoretical studies are 
not conclusive. A probable reason for this discrepancy is that 
most of the theoretical works base their conclusions on crite¬ 
ria which only take into account the brightness distribution 
of the synchrotron emission. 

In the present work, we explore if a theoretical polar- 
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Figure 4. Synthetic polar-referenced angle maps obtained for the 
quasi-parallel (panel a) and quasi-perpendicular (panel b) cases. 
These maps were rotated 60 degrees with the purpose of compar¬ 
ing them to the observations. The axes are in pc and the linear 
color scale is given in degrees. 


ization study can be a useful tool to shed some light on 
this question. We synthesized maps of the polarized emis¬ 
sion that include Stokes Q and U, and the distribution of 
the polar-referenced angle Xr- In our analysis, we considered 
two scenarios: (1) an equatorial belt model where the bright 
rims of SN 1006 are formed from particle acceleration that 
occurs when the shock passes through interstellar magnetic 
fields that are perpendicular to the shock normal (quasi¬ 
perpendicular); (2) a polar cap model where the bright rims 
are formed from acceleration that occurs when the shock 


© 2002 RAS, MNRAS 000, [I]-?? 



















6 E. M. Schneiter 


passes through magnetic fields that are parallel to the shock 
normal (quasi-parallel). 

The resulting polarized intensity maps alone are not 
sufficient to decide which mechanism is the adequate to ex¬ 
plain the radio morphology of SN 1006 since both scenarios, 
quasi-parallel and quasi-perpendicular, yield similar results. 
On the contrary, the synthetic Q maps for the quasi-parallel 
and quasi-perpendicular cases differ much from each other. 
The quasi-perpendicular Q map displays regions with neg¬ 
ative values resembling the bright rims observed in the po¬ 
larized intensity maps. This result is in disagreement with 
the observations, which mainly display positive values (Fig. 
[3] see also Reynolds fe Gilmore|1993 ). In addition, our sim¬ 
ulated Q map for the quasi-parallel case coincides with ob¬ 
servations in that the maxima are spatially coincident with 
the bright arcs of SN 1006, as obtained in our synthetic Q 
map for the quasi-parallel case. Finally, the synthetic polar- 
referenced angle map for the quasi-parallel model is in very 
good agreement with the results reported by [Reynoso et al. 
(2013). 

Our results agree with Bocchino et al. (2011) in that the 
ISM magnetic field must be parallel to the galactic plane, 
tilted with respect to the los , and with a gradient that ex¬ 
plains the assymetry in the emission between the arcs. The 
difference between their result and ours lies in the degree 
of inclination with respect to the los. For the parameter Q 
to be comparable to the observations, we employed an in¬ 
clination of 75°with respect to the los (15°respect to the 
plane of the sky), which is somewhat larger than the best 
fit (38° =b 4°) reported in Bocchino et al. ( 2011| ). 

We calculated the polar-referenced angle (xr) distribu¬ 
tion and found a good agreement with the observations of 
SN 1006 of Reynoso et al. (2013). The maps of observed 
Xr show small scale structure that is absent on our models 
due to the well ordered field imposed. This technique might 
prove helpful for future studies that include the inhomogene¬ 
ity of the media or of the magnetic field structure. 

In summary, our polarization results provide a strong 
support to the quasi parallel model, which is in agreement 
with previous observational ( Reynoso et al.|20l3 ) and theo¬ 
retical works (Bocchino et al. 2011). More importantly, the 
Stokes parameter Q proved to be a powerful tool to de¬ 
termine the particle acceleration mechanism that best ex¬ 
plains the observed morphology and polarization seen in 
SN 1006. We note that in a previous 2D MHD numerical 
simulation, Schneiter et al. (2010) concluded that the quasi¬ 
perpendicular model was the most suitable to explain the 
observations. Our paper shows the importance of using 3D 
simulations, since limiting the analysis to what is possible 
with an axisymmetric setup led to an opposite result. Only 
a 3D model can account for a magnetic field that is tilted 
with respect to the los and can provide a realistic descrip¬ 
tion of the particle acceleration and synchrotron emission of 
an SNR. 
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